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I. INTRODUCTION 

It is well known that higher-derivative gravity has a scalar degree of freedom in general 
[1-3]. In cosmological models of higher- derivative gravity, the scalar mode is expected to 
play an important role [4-6]. On the other hand, some cases are also known that higher 
order terms in curvatures for a gravitational action do not affect cosmological development 
of a scale factor. For example, it is known that terms which consist of contraction of Weyl 
tensors in a gravitational Lagrangian do not change evolutional equations for a scale factor in 
a model with homogeneous and isotropic space. The other special combinations of curvatures 
are known. In the specific dimension, the Euler form as a Lagrangian does not produce the 
dynamics of gravity at all, because the action becomes a topological quantity in such a case. 

The dimensionally continued Euler densities have also been studied [7-13], because of 
their relation to the effective Lagrangian of string theory, and are found to give no scalar 
mode since the second derivative of the metric disappears in the action if we perform in- 
tegration by parts. The absence of scalar modes is interesting for studying black holes in 
the theory, because the scalar modes lead to singularities, in general, which avoid expected 
horizons. 

In recent years, it turned out that there is a special case where a scalar mode disappears in 
higher-derivative gravity. Originally, this fact was found in research of a three-dimensional 
theory of massive gravity [14] and an extended version in four dimensions was proposed 
[15]. The authors of those papers intended to study the renormalizability and unitarity of 
gravitation theory in a maximally-symmetric spacetime. Thus, the absence of a massive 
scalar mode is at least a necessary condition of such theories referred as critical gravity. 1 
Until now, however, only the cases with curvature tensors of a limited number have been 
investigated in higher dimensions [16, 17]. We are interested in the higher order theory of 
gravitation in which a scalar mode does not appear in a general higher dimensions. 

In the present paper, we generalize the structure of the Lagrangian of critical gravity, to 
models with higher order terms in curvature tensors in higher dimensions. We show that such 
extensions can be attained by use of the Lovelock tensors. In order to offer a systematic way 
to construct the required higher order term, we take an explanatory approach by assuming 
the Friedmann-Lemaitre- Robertson- Walker (FLRW) metric. In this approach, the absence of 

1 In our analysis, we do not care for the other critical value for the cosmological constant, etc. 
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second derivatives of the scale factor from the Lagrangian with appropriate total derivatives 
is considered as a necessary condition for disappearing scalar modes. 

It should be noted that the combination defined in D-dimensional spacetime 

R^R^ - _ ^R 2 , (1.1) 

is used in critical gravities in three and four dimensions [14, 15]. This term can be considered 
as a trace of multiplication of the Einstein tensor and a linear combination of the Einstein 
tensor and its trace part: 

(V - \R9v) R" u - 2(rhi) R!r = G ^ (°" u ~ D~i G ^ r ) ' (L2) 

where = R^ v — \Rg^ is the Einstein tensor. Incidentally, R^ u — 2 {d-i) ^ 9^ * s known 
as the Schouten tensor up to the factor (D — 2) _1 . Now, if the FLRW metric is assumed, 
the time-time component of the Einstein tensor does not include the second time- derivative 
of the scale factor. The Schouten tensor appearing the above term is made for its spatial 
component to have no second time- derivative of the scale factor. Therefore, the trace of the 
product of two tensors is linear in the second time- derivative of the scale factor, and if a 
surface term is suitably assigned, the Lagrangian is expressed only with the scale factor and 
its first time-derivative. Thus, additional scalar modes do not appear. From this observation, 
we find that it can be extended by using the Lovelock tensor instead of the Einstein tensor 
when the dimension of spacetime is higher. In the present paper, we do not analyze the 
massive tensor modes in our models. Thus, the genuine criticality as quantum gravity is left 
for future works. 

The FLRW geometry is known to be conformally flat [18], i.e., the Weyl tensor for the 
FLRW cosmological metrics vanishes. An extension of Lovelock gravity for conformally-flat 
geometry was considered by Meissner and Olechowski [19]. They showed that the extension 
is possible for a (R) n term in D dimensions provided n < D, whereas the Lovelock gravity 
has the (R) n term at most 2n < D (where (R) denotes a general curvature tensor). Oliva 
and Ray also construct higher-derivative gravity with the second order equation utilizing the 
Weyl and Riemann tensors [20]. Their Lagrangian involves up to (R) n term with 2n < D 
because of the use of the generalized Kronecker delta as in the case of the Lovelock gravity. 

In the present paper, we show that it is possible to continue the higher-curvature terms 
beyond the number of dimensions. 
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The outline of this paper is as follows. In §2, we construct the candidate Lagrangian 
for cosmological models without scalar modes in the tensorial form as the first step. The 
confirmation of the property of the model Lagrangian is performed in §3, substituting the 
FLRW metric. In §4, we show that an extension to more higher order terms in curvatures 
can be obtained. Using the higher order Lagrangian so far obtained, we propose (toy) models 
for the scale factor with inflationary behavior is shown in §5. In §6, the application of our 
Lagrangian to constructing domain wall solutions is studied. In the last section, we offer 
some concluding remarks. 



II. LOVELOCK TENSORS AND GENERALIZATION OF THE HIGHER 
ORDER TERM IN CRITICAL GRAVITY 



In this section, we construct the higher order term in curvatures by generalizing that of 
critical gravity. We will verify the absence of scalar modes in cosmological models with the 
terms in the next section. First, we introduce the dimensionally continued Euler density 



7" (n) 9 — n z<TiTi---<T n T n rjAipi pA„p n 



(2.1) 



where the generalized Kronecker delta is defined as 



V\ V 2 



AMp AMp 
vi v 2 



St 1 ," 



(2.2) 



The dimensionally continued Euler density consists of n-th order in the curvature tensors 
((R) n ). For example, for n — 1, we find the Einstein-Hilbert term 



(2.3) 



and for n = 2, we find the Gauss-Bonnet term 



L (2) = RavoaR^ - ^RavR^ + ^ , 



(2.4) 



as is well known. 
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Next, we consider the Lovelock tensor [7]. The Lovelock tensor is a generalization of the 
Einstein tensor, and defined as, 

r<(n)V — \ H/j ^\f-~9^ ^) up _ 9 -(n+l) r/xo-m-cr„r„ pAipi pA„p„ /q r\ 

U ^ = 5 ~~ ^ipi-A„p„- rt <Tiri -K a n r n • l^-OJ 



-5 V 

This is a symmetric tensor of ra-th order in the curvature tensors. For example, for n — 1, 
we find 

GS=iV-^%^, (2.6) 
and this is known as the Einstein tensor especially. For n = 2, we find 

= ^{R^tRu PC7T — ^RfipucrR^" 7 — ZR^pRv + RRpu) — -L^g^ u . (2.7) 



It should be noted that 



GW = = - , (2.8) 



where .D denotes the dimension of spacetime. 

Here, we construct the new combination of the Lovelock tensor and the metric multipli- 
cated by the trace of the Lovelock tensor. That is, 

Sft = G^ u - j^rjG^g^ . (2.9) 

It is worth noting that 

GJ = (W»-^)?)' (2-10) 
For example, for n — 1, we obtain 

Sft. = W v - — — ^ -Rg^ , (2.11) 

(i) 2(D-1) 

which is proportional to the Schouten tensor. Therefore we obtain the following combination 

G$S$ = R^BT - J^Tf 2 , (2-12) 

which appears in critical gravities [14, 15]. 

Now, we find that the natural generalization of this is given by 

It is worth pointing out that the expression is symmetric against the exchange of n and n' . 

In the next section, we confirm that this combination is suitable for an extension of 
critical gravity in higher dimensions, by utilizing the FLRW metric. 
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III. HIGHER ORDER TERM FOR FLRW METRIC 



We consider the following FLRW metric in D dimensions: 



ds 2 = -dt 2 + aitydVL^^ , 



(3.1) 



where a(t) is the scale factor and dQjj^ denotes the line element of a maximally symmetric 
space of (D — l)-dimensions, whose scalar curvature is normalized to (D — 1)(D — 2)k with 
k = 1,0,-1. 

First in this section, we examine the Lovelock tensors of (R) n . By explicit calculation of 
curvatures, we find, for n = 1, 

A2 



(1) o (D-l)(D-2) d 2 + k 
G °" 2 a 2 ' 



(D _ 2)- + P- 2 )P-3) « 2 + fe 



(3.2) 
(3.3) 



where a = % and d 



d?a 



. Here and hereafter, we use the suffixes denoting the spatial 



dimensions, i,j — l,2,...,D — l. Also, for n = 2, we obtain 



G (2)0 = p-l)p-2)(D-3)(D-4) / d 2 + ^ 



3 



2(D - 2)(D - 3)(D - 4) 



d d 2 + /c 



a cr 



+ 



(/J - 2)(D - 3)(D - A)(D - 5) fa 2 + k^ 



(3.4) 



(3.5) 



By the combinatorial property of the generalized Kronecker delta, we can find the Love- 
lock tensor for a general n as follows: 

(D - 1){D - 2) •••(£>- 2n) (a 2 + k\ n 



r (n)° 



n(D-2)(D-3)...(£)-2n)- ' " 



n— 1 



a \ cr 



+ 



(£> - 2)(D - 3) • • • (£> - 2n - 1) /d 2 + A; s 



(3.6) 



(3.7) 



and G^ , = G^q = 0. It should be noted that the 00 component of the Lovelock tensor 



does not include d. 
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Next, we calculate the generalized Schouten tensor S( n fi for the FLRW metric. Because 
the trace of the Lovelock tensor is given as 



1 



D — 1 



_ G {r, 



n(D-2)(D-3)---(D-2n 



a ( a 2 + k 



n-l 



a \ a 



+ (D — 2n) 



(D - 2){D - 3) • • • (D - 2n) a 2 + k 



(3.8) 



we find 



/ * 2 1 

S in) ° = n(D - 2)(D - 3) • • • (D - 2n)- ' " ' 



n-l 



a \ cr 



_ {2n _ 1) (D-2)(D-3)---(D-2n) f^d 2 + kj 
% (D-2)(D-3)--(D-2n) ( 'a 2 + A:\" 



(3.9) 
(3.10) 



By construction, a is absent in S(n).. 

Now, we consider the combined term G {jl \ u S( yn >^ 1 ' . We obtain 

G^^f = --AD - 1)[(D - 2){D - 3) ..•(£> - 2n)\[{D - 2){D - 3) • • • (D - 2n')\ 



x 



2(n + n) 



, a ( a 2 + k\ 



n+n'-l 



a \ a 2 J 



+ (D-2(n + n')) 



'a 2 + k 



n+n' 



(3.11) 



Because this combination is apparently linear in a, the action including this term can be 
expressed as the functional of a and a, by means of integrations by part. Therefore, we 
realize that there is no scalar mode in a cosmological setting. 

The combination = — AG^^S^i)^ can exist for limited numbers (n, n'), which 

depends on D. For the exchange symmetry, we assume n < n' . We find that L^ n ^ n '-* is 
non-zero when (n, n') = (1,1) for D — 3,4, (n,n') = (1, 1), (1, 2), (2, 2) for D = 5, 6 and 
(n,n') = (1, 1), (1,2), (1,3), (2,2), (2,3), (3,3) for D = 7,8, etc. Therefore, we have new 
higher-derivative terms including (R) m with m < D — 1 for odd D and m < D — 2 for even 
D. 

Up to now, we become aware of a similarity to the Lovelock Lagrangian. It takes the 
form for the FLRW metric as follows: 



L (m) = (£> _ !)(£> - 2) •••(£>- 2m + 1) 



X 



2m 



a I a + k 



m—l 



a \ cr 



+ (£> - 2m) 



'd 2 + k" 



(3.12) 
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The equivalence up to the overall constant is obvious, that is L oc L^ n+n '\ Incidentally, 
we can consider L(°)( n ) as the Lovelock Lagrangian L^ n \ We note, however, that the Lovelock 
Lagrangian has its meaning as a part of Lagrangian when m = 1 for D — 3, 4, m — 1,2 
for D = 5, 6 and m — 1, 2, 3 for D = 7, 8, etc., because IA m ' becomes a total derivative for 
D = 2m. 

Before closing this section, discussion on relation to the work of Meissner and Olechowski 
[19] is in order. Their approach is equivalent to considering the Lagrangian constructed from 
the Lovelock Lagrangian in which the Riemann tensor is expressed by the Schouten tensor 
under the assumption of vanishing Weyl tensor. 2 They used the combinatorial property of 
the generalized Kronecker delta for extension to higher order in curvatures. If the FLRW 
metric is substituted, we find that the Lagrangian of Meissner and Olechowski of order of 
(R) m coincides with our Lagrangian L( n )( n '), where n + ml = m, up to the constant factor. 
The allowed spacetime dimension is the same for odd m, m < D. Therefore, our Lagrangian 
and theirs are almost equivalent. The variety with respect to two integers in is due 

to the use of Riemann tensors as well as scalar curvatures and Ricci tensors in our approach. 
It is notable that differences may occur if we consider the black hole or non-conformally flat 
solutions in the theory governed by the Lagrangians of higher order terms. 

Later, we show that the restriction by the dimensions can be overcome. To exhibit the 
discussion on the subject, we examine the cosmological action in the present model again in 
the next section. 



IV. FURTHER EXTENSION TO HIGHER ORDER IN CURVATURE 

(ESPECIALLY FOR m > D) 

We consider the action for m > 2: 

/m 
d D x^~g £ a (n)(m _ n) i (n)(m - n) , 

n=0 

with arbitrary coefficients (XuAr m -n)- Here, we regard L^™) as L^ n \ 
Then, the action for the scale factor a(t) can be read as 



S {m) oc f3 m / dta ' 1 



a ( a 2 + k 
2m- — 



m— 1 



+ (D — 2m) 



a 2 + k 



(4.1) 



(4.2) 



a \ a* J \ cr j 

2 Strictly speaking, they replace Riemann tensor by the Kulkarni-Nomizu product of the Schouten tensor 
with the metric in the Lovelock tensor. 
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where 

Pm = (D - l)0!( n )( m _ n ) 

Using the expansion 



(D-2)\ 



(D-2)\ 



m— 1 

(a 2 + k) rn ~ l d = £ 



m — 1 

£ 



(D — 2n— 1)! (D - 2{m - n) - 1)! 

\ m-l um—l—f ( 



k 



m-\-t-2t 



2£+ 1 



V 



i 



dt 



a 2i+1 



(4.3) 



(4.4) 



we rewrite the action, after partial integration, as 

m—l fom—l—f ( 



I (m ^[a,a] = (D — 2m)[3 m / dta^ 1 



-2m ^ 



2£+ 1 



m — l 

I 



a*+* / g2 + £ 



,2m 



(4.5) 



Especially for /c = 0, we find a simple action 

jM = (D-2m)(3 m J dta - 1 



2m — 1 \a? ) 

Now, we define the Lagrangian L^ m ^(a,a) for the scale factor a(t) by 



/M = j dt\S m \a,d). 



(4.6) 



(4.7) 



Then, we find the equation of motion for the scale factor 
d <9L (m) <9L (m) 



dt da da 
= -(D-2m)f3 m a D - 2 



2m— 



d ( a 2 + k 



m—l 



a \ a" 



+ (D-2m- 1) 



a 2 + k 



(4.8) 



The Hamiltonian constraint is regarded as the variation equation = 0, where N is the 



lapse function defined as dt = N(t')dt'. We set N = 1 after the manipulation. We now find 
the Hamiltonian constraint equation 



ON 



(D - 2m)f3 m a 



D-i (a 2 + k s 



(4.9) 



N=l 



It is known that the variation of the lapse function corresponds to the variation of goo 
and the variation of the scale factor corresponds to the variation of ga, up to certain factors. 
We can find, from (4.9) and (4.8), the following generalized Lovelock tensor as the variation 
of the action with respect to the metric: 



G (m)0 = 



)0 _ (D - 2m) (3 m ( a 2 + k s 



(D - 2m)(3 n 
2(D-1) 



2m 



a ( a + k 



m—l 



a \ a A 



+ (D-2m- 1) 



a 2 + k 



(4.10) 
5). (4.11) 
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The generalized Lovelock tensor has a tensorial form, which is proportional to 

1 xq(m) 



(4.12) 



as seen from the construction, but with arbitrary coefficients in the definition of the action 
g(m) j n gpj^g f th e arbitrariness, the functional form of the Lagrangian is unambiguous if 
the conformally flat metric is substituted. 

In the similar manner, the corresponding generalized Schouten tensor is defined as in the 
previous section, 



qflf 
V) 



1 



D 



-G {m) g» u , 



with G^ = g^G {m ^ v . Then we find 
(D - 2m)/3. 



Q _ 
J (m)0 — 



S(m)) 



a a 
2m- - 
a \ 



m—l 



2(D-1) 

(D - 2m) (5 m f a 2 + k 
2(D-1) [ a 2 



a' 

m 



(2m - 1) 



'a 2 



(4.13) 

(4.14) 
(4.15) 



Then, the trace of the product of these tensors is found to be 



L (m)(m') = _ A Q{m)§^ 

(D-2m)(D-2m')f3 m f3 n 
D — 1 



2(m + m')- + {D - 2{m + m')} — 
a 



k 



k 



m+m' — l 

I (4-16.) 



Again, we obtain the same functional form as those constructed from the Lovelock tensors 
in the previous section. More iterative operations yield the more higher order terms. For this 
time, however, no limitation on the relation between the number of dimensions D and the 
order of curvatures I = m + m', except for the case with 21 = D (the Lagrangian becomes a 
total derivative in this case for the conformally flat spacetime). Although we do not exhibit 
an explicit tensorial form of the higher order term in mass dimensions, obviously it can be 
expressed as the combination of curvature tensors and their covariant derivatives. 

Unfortunately, only for D — 4, our iterative method cannot give a (R) 3 term, since 
the (R) 2 term becomes a total derivative for the conformally flat metric. We define the 
expression for by the manner followed by Meissner and Olechowski as 



L {Z) ^^S {1 ZS {1 ^S {1) 



V3 
fJ-3 



(4-17) 

Incidentally, it is obvious that is at most linear in a if the FLRW metric is substituted, 
because only £(1)0 is the component of the Schouten tensor which is linear in a. Then, the 
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desired terms of all the order in curvatures and derivatives in any dimensions can be created 
by our iterative method. 

Another way to cross the dimensional limitation, which is inspired by the work of Meiss- 
ner and Olechowski [19], is also found. In the manner of their paper [19], the following 
Lagrangian (modulo volume form) is proposed: 

where it is expressed by our present notation. We become aware of an extension of it as 

KV^-Sw:, (4-19) 

where mi (i — 1, . . . , n) are arbitrary integers. For the FLRW metric, this term turns out 
to be proportional to 

kB^T^M^)'- (4 ' 2o) 

where £ = YJi=\ m i- This is also the same form of the candidate Lagrangian. Many equivalent 
combinations exist for higher order terms under the vanishing- Weyl-tensor condition. 

In the next section and after, we will turn our attention to apply the higher order La- 
grangian obtained here to cosmological models . 



V. f{H 2 ) COSMOLOGY 



In this section, we investigate the possible inflationary stage in evolution of the universe 
in the model with higher order terms discussed in the previous sections. 
Let us consider the general cosmological action 



D-l 



2£- 

a 



a ( a 2 + k 



+ {D- 2£) 



k " 



(5.1) 



with an appropriate total derivative term with respect to the time, which removes the second 
derivative of the scale factor with respect to time. 

As usual, the energy momentum tensor of matter is taken as 



12 



with the action for matter S matter - The non-zero components of the energy momentum tensor 
of matter is assumed as T ° = — p and T l - — pSj, where p is the energy density and p is the 
pressure, for the FLRW universe. As is seen in the previous sections, the equation of motion 
derived from the action (5.1) coincides with the linear combination of the components of 
the Lovelock tensors in functional form with the assumption of the FLRW metric. Thus, 
the energy conservation 

p+(D-l)-(p + p) = 0, (5.3) 

(X 

holds generally. This fact is due to the absence of the scalar mode and the fact that we did 
not need to rescale the metric. 

We here give a few model in four-dimensional spacetime with focus on the possibility of 
an inflationary growth of the scale factor. Furthermore, we take k = 0, i.e., assume the flat 
space. Then, the action (5.1) is equivalent to 

A-2£ (a 2 \ 



L(a, a) = a 3 J2fie 



2£ 



(5.4) 



The 00 component of the equation of motion reads 



™A-2l (a 2 
1=0 z 



a 2 =P - 



(5.5) 



If we can choose the coefficient freely, almost arbitrary equations including p and 
H = - can be made. That is, 



f{H 2 ) = 



8nG 



(5.6) 



where f(x) is a function which can be expressed by a series and does not include the x 2 
term. Here the factor given in the right hand side is chosen as for the similarity with the 
standard cosmology. We wish to call the cosmology of the model L f(H 2 ) cosmology'. 3 
Now, we specify the function f(x). For example, let us consider 



"{'-Mi) =—"• 

with a typical mass scale M. Solving the equation for H 2 , we obtain 

-l 



(5.7) 



H 2 



8nG 



-P 



1 



+ 



+ 



/8vrG7A 2 
V 3M 2 ) 



(5.8) 



3 If k 7^ 0, we obtain \f{H 2 + k/a 2 ) cosmology'. 
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If the energy density is sufficiently low, such as p <C G 1 M 2 ) the relation in the standard 
cosmology is valid: 

H^^p. (5.9) 
On the other hand, in the era of the high energy density such that p G~ l M 2 , we find 

H 2 ^M 2 , (5.10) 

and we obtain an approximate de Sitter inflation. 

Another model can be chosen, in which the correction is a monomial of the higher order 
term. That is, 

H 2 + M 2 (^Y= 8 -^p (m>3). (5.11) 



M 2 J 3 

In the rapid expansion phase, H 2 3> M 2 , the correction would be dominant. Further, if we 
assume the dust matter, i.e., p oc a -3 , we obtain 

a(t) oc rr . (5.12) 

This indicates the phase of the power-law inflation. 

Even though we show special toy models here, we find that the existence of higher order 
terms yields the inflationary growth of the scale factor with ordinary matter, and with no 
scalar mode and no redefinition of the metric. 

VI. DOMAIN WALL 

We shall apply the previous discussion on a conformally flat metric to solutions for domain 
walls ("thick branes"). The D-dimensional metric suitable for a co-dimension one domain 
wall is given by 

ds 2 = e 2A{y) {-dt 2 + rfx 2 ) + dy 2 , (6.1) 

where ). In addition, we consider a neutral, minimally-coupled, self- 

interacting scalar field as a classical matter field. It is known that superpotential method 
[21-25] is available in this case to find BPS kink equations. Let us investigate the case with 
higher order terms in a similar manner. 

Substituting the metric (6.1), the gravitational Lagrangian becomes 

L = e {D - 1)A f;(-l)'ft [2tA" (A'f~ 2 + D (A'f] , (6.2) 



14 



(A' 



2£ 



where the prime ( ') denotes a derivative with respect to y. The term in A" can be removed 
by discarding a total divergence. Thus we obtain 

L(A,A')^e^ A ±(-lYp e \- D ~ 2i 

The action for the real scalar field <p{y) for static branes is written by 

1 



(6.3) 



MA,0,0') = e (D - 1)A 



(6.4) 



Since the coefficients f3 £ can be arbitrarily chosen, we use an arbitrary function F(x) to 
represent the general action. Then, we can write the total action as 



L toto KA^,0,0') = e (D - 1)A 



F{A l2 )--^'f-V^) 



(6.5) 



The field equations can now be derived in a usual manner. The differential equation for 
the scalar is 

0" + (£-l)Ay-V^ = O, (6.6) 



where — The equation for A is written as 



2[F X (A' 2 ) + 2A' 2 F XX (A ,2 )]A" -(£>-!) 



F(A' 2 ) -2A' 2 F X (A' 2 ) - -(0') 2 -V{<j>) 



0, (6.7) 



where F x = dF J- x ^ and F xx = d . The reparametrization invariance of y leads to the first 

ax ax 1 ,y 

integral given by 

F(A' 2 ) - 2A' 2 F X (A' 2 ) + i(0') 2 - V{<j>) = . (6.8) 

Define h(x) = 2xF x (x) — F(x). Then, one can find h x (x) = jr = F x (x) + 2xF xx (x). 
Note that h(x) as well as F(x) does not possess the x D l 2 term for even D. Using these, the 
equations can be rewritten as 



2h x (A' 2 )A" + (D-l)(<j ) / ) 2 = 0, 



(6.9) 



and 



-/,(A' 2 ) + -(0') 2 -n0) = O. 



(6.10) 



Now, we take a BPS ansatz [21-26] 



A' = -W{4>) 



(6.11) 



15 



Then, the equations reduces to the first order equation 



and the potential 



<!>' = ^Yh x (W 2 )W^ (6.12) 



V(<P) = p^[/^ 2 )] 2 (Wg 2 - h(W 2 ) , (6.13) 



in terms of W{4>). 

If the function h(x) is a polynomial up to the quadratic order, a domain wall solution 
exists, for the potential V(4>) has two minima [21-26]. The solution can be expressed as a 
kink solution in the 4 theory. We find here that the potential with many vacua naturally 
corresponds to the general higher order gravity. 

For example, we try to express the sine-Gordon equation. We take W(4>) = B<fi, with a 
positive constant B. Further, we choose 

h x (x) = Ccos^-, (6.14) 

with constants a and C. This choice is possible if the spacetime dimension D is odd. Then 
the scalar obeys the sine-Gordon equation 

aB 2 C aB<j> aB 3 C 2 nl . 

(f) = ~ £)_ 1 sm ~2~'P = ~ ( D _iy amaB< l > > ( 6 - 15 ) 

and an exact static solution is known as 

(f)(y) = ^arctane^^^ — ^ . (6.16) 

Then, the potential takes the form: 

\r/j.\ 2£ 2 C 2 2 aB<P 8C (aB<P . aB<P ^B^ 
V((p) = — — - cos — — — sin — h cos — — + V . (6.17) 



(D-l) 2 2 a 

The minima of the potential are located at aB<p/2 = tt/2 ± nn (n = 0, 1, 2, . . .). Although 
it is difficult to obtain exact solutions of other types, we can suppose multiple domain walls 
with distinct topological numbers in this and similar models. The model with many vacua 
may also serve an interesting mechanism to realize naturally a small cosmological constant 
in (thick) brane world. This possibility will be studied in future. 
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VII. SUMMARY AND CONCLUSION 

In the present paper, we have attempted to show the possible quasi-linear second-order 
theory of gravity in conformally flat spacetimes. Models with arbitrary higher order of cur- 
vatures have been obtained. As long as we adopt an isotropic and homogeneous cosmological 
setting, the energy conservation holds in the models because there is no scalar mode and no 
requirement of frame rescaling. In spite of them, inflationary expansion can be found in the 
models. We have also found that the domain wall solution in the present type of the higher 
order gravity can be obtained naturally with the potential having many minima. 

Our work corresponds to the extension of the Lovelock higher-curvature gravity in arbi- 
trarily higher order terms in higher dimensions. Our analysis, however, has been limited for 
the case with the conformally flat spacetime and the coefficient on the tensorial form of the 
action has been still ambiguous. The stability of the solutions obtained here is problematic 
for anisotropic perturbations or tensor modes. To study it, we should classify the tensorial 
form of the Lagrangian which is equvalent only if the Weyl tensor vanishes. It is known that 
the dimensionally continued Euler forms have the property of factorization in terms of those 
of lower orders when the spacetime is represented by the direct product of spaces [10]. Thus, 
compactification should be worth studying to seek some special combination of curvature 
tensors. It is expected that the hopeful 'critical' relation among the coefficients of different 
orders of curvatures will be selected by consideration on various background spacetimes. 

Nevertheless, we emphasize that our arbitrarilty high order gravity can be applied to 
many models in various contexts. The problem of initial singularity can be reconsidered by 
studying classical bouncing universes in our model. On the other hand, the Wheeler-DeWitt 
equation of the Lagrangian should lead to higher-derivative quantum cosmology. The equa- 
tion must be difficult to treat with, but the study on it may shed new light on quantum 
gravity. The possible black hole solutions are interesting in both cases of asymptotically flat 
and asymptotically AdS spacetime. We shall return to some of the problems in future. 
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